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An updated global analysis within the Standard Model (SM) of all relevant electroweak precision 
and Higgs boson search data is presented with special emphasis on the implications for the Higgs 
boson mass, Mr- Included are, in particular, the most recent results on the top quark and W boson 
masses, updated and significantly shifted constraints on the strong coupling constant, a 3 , from r 
decays and other low energy measurements such as from atomic parity violation and neutrino deep 
inelastic scattering. The latest results from searches for Higgs production and decay at the Tevatron 
are incorporated together with the older constraints from LEP 2. I find a trimodal probability 
distribution for Mh with a fairly narrow preferred 90% CL window, 115 GeV < Mh < 148 GeV. 
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I. INTRODUCTION 

One of the prime missions of the Large Hadron Col- 
lider (LHC) at CERN is the search for the Higgs boson. 
Within the SM, its existence is solidly predicted but only 
some semi-quantitative theoretical constraints exist for 
its mass. If the SM is the correct low-energy theory only 
up to a new physics scale which is itself not much larger 
than Mh, one would find at one-loop order and neglect- 
ing all couplings other than the Higgs self-coupling A, the 
" triviality" condition, 
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M H < 816 GeV, 



where v — 246.22 GeV is the Higgs vacuum expectation 
value. This is very close to the bound obtained from the 
requirement of unitarity of the partial S-wave amplitude 
of elastic Goldstone boson scattering [T], 
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M H < 781 GeV. 



Requiring the absence of a Landau pole in A up to the 
reduced Planck scale, up = 2.4 x 10 18 GeV, yields in the 
same approximation as above, 
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M H < 147 GeV, 



while a refined analysis including top quark mass, m t , 
and two- loop effects gives Mh ^ 0(m t ) [2]. Vacuum 
stability, i.e., the requirement that the scalar potential is 
bounded below, implies the lower bound (driven by the 
large top quark Yukawa coupling), M H ^ 130 GeV [3J, 
but the vacuum would still be sufficiently long-lived 
(metastable) for M H > 115 GeV @]. 
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Quantum loop corrections affecting the multitude of 
electroweak precision observables — most importantly 
the W boson mass, Mw, and effective weak mixing an- 
gles, 6yy — allow to simultaneously test the model and 
over-constrain its free parameters including Mh- More- 
over, by comparing various cross-section measurements 
with the SM prediction as a function of Mh one may 
identify preferred and disfavored regions. In this way, 
CDF and D0 at the Tevatron [5] concluded that the win- 
dow, 162 GeV < Mh < 166 GeV, is incompatible with 
their data at the 95% CL. Likewise, the LEP 2 Collabora- 
tions [6] found the 95% CL lower limit, M H > 114.4 GeV. 

In this communication, I update the global analysis of 
Ref. [7] where the statistical method is described in de- 
tail. I incorporate all direct (search) and indirect (preci- 
sion) data, including new radiative corrections and signif- 
icant improvements and changes in several precision ob- 
servables. For alternative analyses, see Refs. [Hill], where 
the latter differs mostly by the neglect of the low energy 
data and the treatment of LEP 2 (see Section |HB[ ). 



II. EXPERIMENTAL DATA 

A. Electroweak precision observables 

The strongest Mh constraints come from asymmetry 
measurements determining sin 2 0$ for leptons [10] at 
LEP 1 and the SLC Q31 Q2], and from M w at LEP 2 
and the Tevatron [T3J . The two most precise determina- 
tions of sin 2 9$ deviate by about 3 a from each other, but 
since both are statistics dominated we consider this in the 
following as a fluctuation. These constraints are strongly 
correlated (39%) with nil, giving great importance to the 
direct determination of the latter at the Tevatron 1141. 



m t 



173.1 ± 0.6 (stat.) ± 1.1 (syst.) ± 0.5 (QCD) GeV. 



This is converted to the MS-mass definition using the 
three-loop formula [15] which gives rise to the QCD er- 
ror (the size of the three- loop term). It is assumed that 
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FIG. 1. la (39.35% CL) contour lines for Mh as a function 
of m t for various inputs and the solid (red) 90% CL ellipse 
(Ax 2 = 4.605) allowed by all data. a s (M z ) = 0.1183 is as- 
sumed except for fits including the Z lineshape or low energy 
data. The lower limit from LEP 2 and the excluded window 
from the Tevatron (both at the 95% CL) are also shown. 



this accounts for the error from interpreting the mass 
extracted at the Tevatron as the pole mass. 

Other Z pole constraints are the Z width, Tz, the total 
hadronic peak cross section, <7had, and a set of branching 
ratios, Ri [TU]. Finally, there is a wide range of low en- 
ergy experiments from atomic parity violation (APV) to 
neutrino and polarized electron scattering. See Ref. [16] 
for tables of inputs almost identically to those used here. 

All experimental and theoretical uncertainties and cor- 
relations are included in the fits. The error from unknown 
higher order electroweak corrections is implemented by 
allowing the so-called oblique T-parameter [TT] to float 
subject to the constraint T = 0±0.02. Errors from differ- 
ent sources have been added in quadrature and in most 
(but not all) cases been treated as Gaussian. The result- 
ing constraints are depicted in Figure [I] Some comments 
on those inputs which have shifted recently or which have 
been discussed controversially are in order: 

The theoretical predictions for Mw and sin 0<§ need 
the renormalization group evolution of the electromag- 
netic coupling from the Thomson limit to the weak scale. 
Entering the implementation (in the FORTRAN pack- 
age GAPP (18j) is the MS definition, a(M z ), which 
is updated from Ref. |19j with its central value moved 
upwards and its uncertainty almost halved. The cor- 
responding hadronic vacuum polarization effects can be 
translated from cross-section data for e + e~ — ¥ hadrons, 
which in turn can be obtained by standard e + e~ annihi- 
lation or by the high statistics (but systematics dom- 
inated) method |20] of using radiative returns from a 
IS resonance. In addition, there are measurements of 
r decay spectral functions which can be included with 
the appropriate isospin corrections |21j . However, the 



results reveal some discrepancies. The r data imply 
lower values for the extracted Mh of about 6% com- 
pared to the e + e~ data. This conflict is smaller than 
in the past and some of it appears to be experimen- 
tal. The dominant e + e~ — > tt + tt~ cross-section has been 
measured by CMD-2 [22] and SND [23] and the results 
are in good agreement with each other, but are lower 
than those obtained from T(AS) radiative returns by 
BABAR [23]. In turn, the latter agrees quite well with 
the r analysis including the energy dependence (shape). 
In contrast, the shape and smaller overall cross-section 
from 7r + 7r~ pairs radiatively returned from the <E> and 
detected by KLOE [23 differ significantly from BABAR 
(a recent review on the e + e~~ data is Ref. [26). All 
measurements including older data are accounted for on 
the basis of results from Refs. [22 [2S1 [21] • The correla- 
tion with the magnetic moment and the non-linear 
a s dependence of a(Mz) are addressed. The correla- 
tion of a(Mz) with a s has been treated by using as 

input (fit constraint) instead of Aa^ d (Mz) the analo- 
gous low-energy contribution by the three light quarks, 

Aa had( L8 GeV ) = ( 57 - 29 ± °- 90 ) x 1(r4 > and b y ec- 
lating the perturbative and heavy quark contributions to 

a(Mz) in each call of the fits according to Ref. [TH]. The 
error is from e + e~ data below 1.8 GeV and r decay data, 
from uncertainties in the isospin breaking effects (af- 
fecting the interpretation of the r data), from unknown 
higher order perturbative and non-perturbative QCD ef- 
fecs; and from the charm and bottom quark masses. 

There is extra information on sin 2 6yy and Mh in the 
Z boson vector couplings, which is used best if a s is 
constrained independently. For this I use the extraction 
of a s from the r lifetime, t t , because (i) the r scale is 
low, so that upon extrapolation to the Z scale, the a s 
error shrinks by an order of magnitude; (ii) this scale 
is still high enough that the operator product expansion 
(OPE) can be applied; (iii) t t is fully inclusive and thus 
free of hadronization effects; (iv) OPE breaking effects 
occur only where they are kinematically suppressed; (v) 
non-perturbative effects can be constrained by experi- 
mental data; (vi) the complete four-loop order (massless) 
QCD expression is known; and (vii) large effects associ- 
ated with the QCD /3-function can be re-summed [2H] m 
contour improved perturbation theory (CIPT). However, 
while CIPT shows faster convergence in the lower calcu- 
lated orders, doubts have been cast on the method by the 
observation that at least in a specific model [2H] including 
theoretical constraints on the large-order behavior, ordi- 
nary fixed order perturbation theory (FOPT) may nev- 
ertheless give a better approximation. The largest uncer- 
tainty arises from the truncation of the FOPT series and 
is taken as the a 4 s term. I find a s (M z ) [r] = 0.1174^ 0°^ 
which updates Ref. [3U]. The effects of using FOPT in- 
stead of CIPT, of using the theoretically better motivated 
spectral functions of Ref. [3T] in place of previous results, 
and of including the four- loop result [52] i all significantly 
reduce the extraced a s value. 

There are precise APV experiments in Cs [53] [M] and 
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FIG. 2. \ 2 distribution of Mh from all data. The horizontal 
solid (green) line marks the \ 2 minimum, Xmin — 39.88, while 
the dashed (blue) lines refer to integer values of yj x 2 ~ Xmin- 
The combined effects of precision, LEP 2, and Tevatron data 
result in the pronounced dip at Mh = 115.8 GeV. Values of 
M H > 141 (158) GeV are excluded at the 2 (3) a level. 



Tl [351 HI] i where the error associated with atomic wave 
functions is quite small for Cs [37]. The extracted weak 
mixing angle (in the MS-scheme), = 0.2314 ± 0.0014, 
now agrees perfectly with = 0.23116 ± 0.00013 from 
the SM fit, where the theoretical effects in Refs. [3"T1 [3"5] 
together with an update of the SM calculation [39] re- 
moved an earlier 2.3 a deviation from the SM. 

Neutrino-nucleus deep inelastic scattering (^-DIS) is 
dominated by the NuTeV result [JD] for the on-shell 
weak mixing angle, Syy = 0.2277 ± 0.0016, which ini- 
tially was 3.0 a higher than the SM prediction, s\y = 
0.22292 ± 0.00028. Since then a number of experimen- 
tal and theoretical developments shifted the extracted 
Syy, most of them reducing the discrepancy: (i) NuTeV 
also measured [41] the difference between the strange 
and antistrange quark momentum distributions, S~ — 
0.00196 ± 0.00143. The effect of S~ ^ on the NuTeV 
value for has been studied in Ref. [12], and the S~ 
above shifts by -0.0014 ± 0.0010. In view of theo- 
retical arguments favoring a much smaller or negligible 
effect, I take half of the above shift as an estimate of 
both the S~ effect and the associated error, (ii) The 
measured branching ratio for K e3 decays enters crucially 
in the determination of the v e (v e ) contamination of the 
v^iyyu) beam. Since the time of Ref. [30] it has changed 
by more than 4 c, so that a move of by +0.0016 
is applied and the corresponding v e {v e ) error decreased 
by a factor of 2/3. (iii) Parton distribution functions 
(PDFs) seem to violate isospin symmetry at levels much 
stronger than expected. While isospin symmetry violat- 
ing PDFs are currently not well constrained phenomeno- 
logically, the leading contribution from quark mass dif- 
ferences turns out to be largely model-independent [43] 
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FIG. 3. Probability distribution of Mh subject to all data. 
The nominal 95% CL exclusion ranges [5] [6] from LEP 2 and 
the Tevatron are also indicated. 

and a shift, 5s 2 w = -0.0015 ± 0.0003 0HS], is applied, 
(iv) QED splitting effects also violate isospin symmetry, 
shifting s^, by —0.0011 3B] with a 100% assigned error 
(the sign is model-independent), (v) The isovector EMC 
effect [47] model-independently reduces the discrepancy, 
shifting s 2 w by -0.0019 ± 0.0006 [35]. (vi) The extracted 
s^y may also shift significantly when analyzed using the 
most recent QCD [35], QED and electroweak [39. radia- 
tive corrections, but their precise impact will be revealed 
only after the NuTeV data have been re-analyzed with a 
new set of PDFs. Remaining radiative corrections have 
been estimated 49 to induce an error of ±0.0005 in s^. 
With these corrections, the global i^-DIS average in terms 
of the weak mixing angle is = 0.2254 ± 0.0018. 



B. Collider searches 

At LEP 2 with energies up to ^/s w 209 GeV, the 
Higgs boson was searched for in the dominant (ss 74%) 
bb decay channel, produced in the Higgsstrahlung pro- 
cess, e + e~ — > ZH . In addition, the H — > t + t~ channel 
7%) was studied for the Z boson decaying into two 
jets. The combination [5] of the four experiments, all 
channels and all i/i values, resulted in the nominal lower 
bound, Mh > 114.4 GeV. However, the combined data 
are neither particularly compatible with the hypothesis 
Mh = 115 GeV (15% CL), nor with background only 
(9% CL). The reason is that the results by ALEPH are by 
themselves in very good agreement with Mh ~ 114 GeV 
(due to an excess in the 4-jet channel) thereby strongly 
rejecting the background only hypothesis, while the re- 
sults based on the other channels and experiments (espe- 
cially DELPHI) are incompatible with any signal. Over- 
all, a signal for 115 GeV < M H < 119.5 GeV is favored 
by the data, but not with high significance. 
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TABLE I. Upper and lower bounds for Mh [in GeV] at various confidence levels (CLs). Shown are the Bayesian limits 
defined as the central intervals of the integrated probability density function (see Figure |3J. Also show (in parentheses and for 
comparison only) are the limits derived from simple frequentist (maximum likelihood) reasoning (see also Figure [2|. The two 
values in the first line are for the median and the best fit, respectively. They differ because the distribution is not symmetric. 
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123.6 


(115.8) 


123.6 


(115.8) 


median (Xmin) 
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137.6 
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95 


5 
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148 




90% CL range 


97.725 


2.275 


114.8 


(114.3) 


155 


(141) 


2cr range 


99.5 


0.5 


114.3 




197 




99% CL range 


99.865 


0.135 


113.9 


(113.5) 


217 


(158) 


3cr range 


99.95 


0.05 


113.6 




232 




99.9% CL range 


99.9968 


0.0032 


112.9 


(112.6) 


276 


(240) 


4cr range 


99.9995 


0.0005 


112.4 




304 




99.999% CL range 


99.999971 


0.000029 


111.5 


(111.3) 


349 


(315) 


5cr range 



The LEP 2 results can be included by adding the solid 
line for the observed log- likelihood ratio (LLR b s ) shown 
in Figure 1 of Ref. [5] to the x 2 -function derived from 
the precision observables in Section II A| The quantity 
LLR bs is defined as —2\uQ(Mh), where Q(Mh) is the 
ratio of the likelihood for the signal of a particular Mh 
hypothesis plus the background to that of the background 
alone. This treatment is rigorous in the limit of large data 
samples and serves as a good approximation otherwise. 
It is emphasized that treating the LEP 2 results as a step 
function with threshold at the nominal lower Mh bound 
is a poor approximation whenever there is a noticeable 
upward fluctuation in the data beyond that threshold and 
results in systematic and significant upwards shifts of the 
upper bounds (compare, e.g., with Ref. [§])• 

At the Tevatron running at y/s = 1.96 TeV, the Higgs 
boson can conceivably be produced in association with 
W or Z bosons, pp — > W/ZH (the counterpart of Hig- 
gsstrahlung at LEP 2), or through gluon (gg — > H) or 
vector boson (qq — > q'q'H) fusion. The studied de- 
cay channels besides H — > bb and H — > t + t~ are 
the one-loop decay H — > 77 as well as (dominant for 
M H > 140 GeV) H -> W+W~ . For their combina- 
tion [50] CDF and D0 analyzed 90 individual processes. 
As was the case at LEP 2, but somewhat more sig- 
nificantly, the low mass Higgs region is favored by the 
data. This is pronounced around 115 GeV, but persists 
until Mh = 155 GeV. On the other hand, the range 
155 GeV < Mjj < 197 GeV is disfavored especially in 
the nominal exclusion window. The Tevatron results are 
incorporated by adding the LLR Q b s column given in Ta- 
ble XIX of Ref. [50] as another contribution to % 2 . 



III. RESULTS 

The following results are based on a Bayesian treat- 
ment [7] which is particularly adequate for parameter 



estimation (as opposed to hypothesis testing). The no- 
torious objection of the necessity of a prior distribution 
notwithstanding, Bayesian data analysis provides a first 
principles setup with strong emphasis on the entire pos- 
terior density [51], It is given by, 



p(M H ) 



LEP 2 ^Tevatron 



M 



H 



(1) 



where the first factor is from the precision data, while the 
factors of Q — Q(Mh) are as described in Section II B 



The last factor is the (improper) non-informative prior 
density chosen such that the variable In Mh has a flat 
prior which one can argue is the most conservative (least 
informative) for a variable defined over the real numbers. 
Alternatively choosing a flat prior in Mh itself (i.e., drop- 
ping the factor M^ 1 ) increases, for example, the 95% CL 
upper limit by modest 3 GeV because small Mh values 
are penalized a priori. The 90% preferred range (95% 
CL lower and upper bounds) for Mh is given by, 



115 GeV <M H < 148 GeV, 



(2) 



with corresponding bounds for other CLs shown in Ta- 
ble [T] The Table also shows that a Higgs boson discovery 
with Mh > 350 GeV would simultaneously mean the in- 
direct discovery of new weak scale physics which would 
have to bridge the gap between the physical Higgs mass 
and the mass derived from the current data when assum- 
ing the validity of the SM. 

Figure [2] shows the distribution of \ 2 = — 2lnp(Mn), 
while Figure [3] represents p(Mn) binned in 1 GeV steps. 
Only by virtue of xfndirect ^ s ^ possible to obtain a 
proper p(Mn) density. A trimodal distribution emerges 
with a tall peak since the searches of both, LEP 2 and 
the Tevatron, see some excess events hinting at M H < 
120 GeV, and also since the high energy precision data 
favor M H = 87±1| GeV [B|. This M H value agrees well 
with Ref. [9, (M H = 80±^ GeV) and the global fit here 
(M H = 91^24 GeV). The other two modes are due to the 
precision data in regions where QTovatron is basically flat. 
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